Spontaneous generation of a half-quantum vortex in spin-orbit coupled Bose-Einstein 

condensates 
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According to the "no-node" theorem, many-body ground state wavefunctions of conventional 
Bose-Einstein condensation (BEC) are positive-definite, thus time-reversal symmetry cannot be 
spontaneously broken. We find that multi-component bosons with spin-orbit coupling provide an 
unconventional type of BEC beyond this paradigm. In a harmonic confining trap, the condensate 
spontaneously generates a half-quantum vortex combined with the skyrmion type of spin texture. 
This phenomenon can be realized for both the cold atom systems with artificial spin-orbit coupling 
generated from atom-laser interaction, and for exciton condensates in semi-conductor systems. 
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The conventional many-body ground state wavefunc- 
tions of bosons satisfy the celebrated "no-node" theorem 
in the absence of rotation, as pointed out by Feynman 
which means that they are positive-definite in the 
coordinate representation. This theorem implies that 
time-reversal (TR) symmetry cannot be spontaneously 
broken. It applies to various quantum ground states 
of bosons including Bose-Einstein condensates (BEC), 
Mott-insulating states, density-wave states, and super- 
solid states, thus making this a very general statement. 

It would be interesting to search for novel types of 
quantum states of bosons beyond this paradigm, such 
as unconventional BECs with complex-valued wavefunc- 
tions. The "no-node" theorem is only a ground state 
property, hence it does not apply to excited states. Re- 
cently, a meta-stable state of bosons has been realized, 
using cold atoms in optical lattices, by pumping atoms 
into high orbital bands (e.g. p-bands) with a long life 
time UQ. Interactions among orbital bosons are char- 
acterized by an orbital Hund's rule 0, H|, which gives 
rise to a class of complex supernuid states by develop- 
ing onsite orbital angular momentum, thus breaking TR 
symmetry spontaneously 0-0]. 

The "no-node" theorem does not apply to spinful 
bosons with spin-orbit (SO) coupling either. The linear 
dependence on momentum introduced by the SO cou- 
pling invalidates Feynman's proof. Although alkali atoms 
are too heavy to exhibit rclativistic SO coupling, artificial 
SO coupling can be induced by laser-atom interactions 
The BECs in such systems have been studied, 
and have been shown to exhibit frustrations of superflu- 
idit y |1 4lj . On the other hand, excitons in semiconductors 
[15l - il8j with small effective mass can exhibit SO coupling 
in the center-of-mass motion [19 21]. In particular, ex- 
citing progress has been made in indirect exciton systems 
in coupled quantum wells where electrons and holes are 
spatially separate 22-25|. The extraordinarily long life 
time of indirect excitons provides a wonderful opportu- 
nity to investigate the exotic state of matter of exciton 
condensation llal. 



In this article, we show that bosons with Rashba SO 
coupling provide a new opportunity to investigate un- 
conventional BEC beyond the "no-node" theorem. We 
find that Rashba SO coupling leads to fragmented con- 
densation in homogeneous systems. Furthermore, with 
an external harmonic potential trap, TR symmetry is 
spontaneously broken. In this case, the ground state 
condensate develops orbital angular momentum, which 
can be viewed as a half-quantum vortex. Moreover, the 
spin density distribution exhibits a cylindrically symmet- 
ric spiral pattern as skyrmions. The experimental real- 
ization, in both cold atom and semiconductor exciton 
systems, is discussed. 

We begin with a 3D two-component boson system with 
Rashba SO coupling in the xy-plane and with contact 
interaction, described by 



H = J d 3 r V4l 



ft 2 V 2 
" 2M 



(1) 



where ip a is the boson operator; the pseudospin indices 
a =~f, J, refer to two different internal components of 
the bosons; \r is the SO coupling strength; g describes 
the s-wave scattering interaction. Although Eq. (JTJ is 
for bosons, it satisfies a suitably defined TR symmetry 
T = 102C satisfying T 2 = — 1, where C is the complex 
conjugate operation and o~2 operates on the boson pseu- 
dospin degree of freedom. 

In the homogeneous system, the single particle states 
are the helicity eigenstates of a ■ (k x z) which have a 
dispersion relation given by e±(k) — jm((^ ~ ^o) 2 + k" 2 ), 
where k = The energy minima are located on 

the lower branch along a ring with radius kg in the 
plane of k z = 0. The corresponding two-component 
wavefunction ijj+(k) with |fc| = fco c & n be solved as 
V>+(fc) = ^(e-^ fc/2 ,ie^ fc/2 ), where 4> k is the azimuth 

angle of the projection of k in the xj/-plane. The in- 
teraction part in Eq. ([T]) in the helicity basis can be 
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FIG. 1: The degeneracy of the single particle ground states 
along a circle with the radius ko in momentum space. We 
assume that the fragmented condensation occurs at points A 
(koe x ,0,0) and B (— ko& x ,0, 0) with spin polarizations along 
Hy, respectively. The low energy excitations within \k—ko | < 
A, \k z \ < A, and A/ko -C 1 are classified into two regimes I 
(inside two cylinders centering around A and B with radius 
of A) and II (outside). 

represented as 

Hint = | X! X! + ^ x \pup)(p-2 -t,iAp2\v) 
\nvpp1p2q 

x VaCpi + Q)i>t(p2 - q)^{P2)^ P {pi), (2) 

where the Greek indices A, v, p, p are the helicity indices 
±; the matrix elements denote the inner product of spin 
wavefunctions of two helicity eigenstates at different mo- 
menta, e.g., (pi + q; A|pi; p) = |[1 + A/5e j(< ^>i~ , ^>i+'i ) ]. 

At the Hartree-Fock level, the kinetic energy can be 
minimized through coherent or fragmented condensation 
in the linear superpositions of the single particle ground 
states. For the interaction energy, the exchange part 
among bosons is positive within the s-wave scattering ap- 
proximation where g is independent of momentum, thus 
the ground state favors the zero exchange energy. The 
spin polarizations at k and — k in the plane of k z = are 
orthogonal to each other. Thus, without loss of gener- 
ality, the fragmented condensate defined below satisfies 
this requirement 

$0 = J— WUk A )] N *Wi(k B )] NB \0), (3) 

where Jca = (— &o, 0, 0) and ks = (ko, 0, 0); and (Na, Nb) 
is the particle number partition satisfying Na+Nb = Nq, 
with Nq the total particle number in the condensate. In 
homogeneous systems, fragmented condensates with dif- 
ferent partitions do not mix because of momentum con- 
servation. 

An "order from disorder" calculation on the low energy 
excitation spectra is needed to determine the optimal 
particle number partition. The momentum cutoff HA is 
defined as h 2 A 2 / (2M) as gn, where gn is the interaction 
energy scale and n is the particle density. We only con- 
sider the limit of strong SO coupling, ko ^> A, and leave 



the general case for later study. The low energy Bogoli- 
ubov excitations that significantly mix particle and hole 

operators mainly lie in the regime | ^Jk 2 + k 2 — ko | < A 

and \k z \ < A, which can be further divided into two parts 

I and II as depicted in Fig[T] Part I is inside two cylin- 
ders with the radius of A centering around A and B, and 
part II is outside these two cylinders. 

For the low energy states in part I, we define 
the bosonic operators in the lower branch as — 
ip + (-k e x +q) and b$= ip + (k e x + q), and (a q = ) = ^JW2 
and (bq—o) — VNb respectively. The relative phase be- 
tween two condensates can be set to zero due to a suit- 
able translation operation. The low energy excitations in 
this region have been calculated in Ref. [14]. By defin- 
ing j\(q) = ^(V^oJ + VN~ b b\), 4{q) = ^(^oj - 
y/Nabjj), the mean- field Hamiltonian in region I, up to 
the quadratic level of q, can be represented as 

Hmfa = £>($)W($) 7 ($)+72($) 72 ($)} 

+ 9N {4(q)ll(-q)+h.c.}, (4) 

where E(±q) = h(q 2 x + q 2 z )/(2M). The 71 mode af- 
ter the Bogoliubov diagonalization yields the spectra 

u% = V / W+2#o) « V^^F+e This is 
the phonon mode describing the overall density fluctua- 
tions, which exhibits linear dispersion relation for (fin the 
xz-plane and becomes soft for q along e y . The 72 mode 
represents the relative density fluctuations between two 
condensates which can also be considered as spin density 
fluctuations. This mode remains a free particle spectrum 
E(q). Both of the 71,2 modes, at q <C ko, only depend on 
the total condensation density Na + N B = Nq. Hence, 
the contribution from part I does not lift the degeneracy 
between different partitions of (Na, Nb) to the quadratic 
order of q. 

Now, we turn to the low energy excitations in part 

II where the dispersion of ipj: , is degenerate with the 
state oi ib r but not with ib„r r , • The mean field 
Hamiltonian reads 

H M f,2 = 5^4 i+ Vjg i+ {e(A;) + |(JVb-AJVcoB0 fc )} 
fe 

+ .9(4 + ^U+VNaN b cos^e^s + h.c.}, (5) 

where A A" = Na — N B - The Bogoliubov spectra can be 
solved as H MF . 2 = Y^k { w (^)(t](^)73(^) + 5)} witn 




After summing over k around the ring, the second term 
in Eq. [5] disappears. Thus the ground state zero point 
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FIG. 2: A) The radial density distribution of spin up and 
down components, and the total density distribution in the 
unit of No at a = 4 and /3 — 40 in a harmonic trap. B) The 
spin density distribution along the a>axis which spirals in the 
z-x plane at an approximate wavevector of 2ko, whose value 
at the origin is normalized to 1. The spin density in the whole 
plane exhibits the skyrmion configuration. 



motion from Eq. [6] selects the optimal partition of AN 
0, i.e., an equal partition between two condensates. 

Fragmentation of BECs is usually unstable against 
perturbations [ 



27] 



It can be removed by non- 
conservation of momentum induced by spatial inhomo- 
geneity. Now we consider a realistic system with confin- 
ing trap. Our main interest is the interplay between the 
confining potential and the Rashba SO coupling which 
only lies in the a;y-plane. Thus we neglect the effect of 
the trapping potential along the 2-axis, and assume that 
the condensate wavefunction is uniform along this axis. 
We consider V ex {r) — ^Muj\{x 2 + y 2 ), and define the 
characteristic SO energy scale for the harmonic trap as 
E so = hXftJl where I — y/h/{MujT) is the length scale of 
the trap; correspondingly, we also define the dimension- 
less parameter a = E so /(FujJt) — IJcq. Due to the 2D ro- 
tational symmetry, the single particle wavefunctions can 
be denoted by the total angular momentum j z — rn + | . 
The ground state single particle wavefunctions form a 
Kramer doublet as represented in polar coordinates as 



4>i(r, <j),z) 



f(r) 



ijj_i{r,<p,z) = Ttpi. (7) 



where both f(r) and g(r) are real functions. 

Let us gain some intuition about f(r) and g(r) by con- 
sidering the strong SO limit a 3> 1. The Fourier com- 
ponents of ip±± in momentum space mainly distribute 
around the low energy ring. They can be approximately 



represented as 



dcj) k / dk w{k) e ±l ^ k \%l; + {k)), 



(8) 



where the weight w(k) mainly distributes along the low 
energy ring with a broadening at the order of 1//. In 
other words, the harmonic potential, which reads V ex = 
^Mwy.(ic?£ — A(k)) 2 in momentum representation, where 
A(k) = i(tp + (k)\d^\ijj + (k)) . For the Rashba coupling, the 
field strength of A(k) corresponds to a 7r-flux located at 
k = (0,0). V ex breaks the degeneracy of the plane wave 
states with \k\ = ko by quantizing the orbital motion 
along the ring as 



1 , , 2 ,m + 

MuU—r- 
2 fc 



i 

2 \2 



_L^ t(to+ I)2. (9) 



The ground state forms the Kramer doublets correspond- 
ing to m + | = ±|, and the low energy azimuthal ex- 
citations are with m+i = ±|,±|.... The excitations 
in the radial direction of change the node structures 
of f(r) and g(r), which cost the energy at the order of 
hwT and thus are much higher energy states than the 
azimuthal excitations at a 3> 1. For the strong SO cou- 
pling a 3> 1, the single particle ground doublet states 
have nearly equal weight in the spin up and down com- 
ponents, i.e. / drdcf) r\f(r)\ 2 rj J drdcf> r\g(r)\ 2 , thus the 
spin moment averages to zero. When the bosons con- 
dense into one of the TR doublets, the average orbital an- 
gular momentum per particle is ft/2, i.e., one spin compo- 
nent stays in the s-state and the other one in the p-state. 
This is a half-quantum vortex configuration 28l430j , thus 
spontaneously breaking TR symmetry. 

Next we turn on interactions and define the character- 
istic interaction energy scale Ei nt = gNo/(irl 2 ) and the 
dimensionless parameter /3 = Ei nt / (Hu>t)- We numeri- 
cally solve the Gross-Pitaevskii equation 

r ft 2 v 2 

I 2M~ + K-Wv°x,ap + iVxOV.a/s) + sC^V^) 

+ ^Mu: 2 r 2 }Mr^)=EMr,4>) (10) 

in the harmonic trap with parameters a — 4 and j3 — 40. 
The interaction effectively weakens the harmonic poten- 
tial and does not change the orbital partial wave struc- 
ture of the wavefunctions. We show the radial density 
profiles of both spin components l/MI 2 and |g(r)| 2 in 
Fig. [2] A. Each of them exhibits oscillations at a pitch 
value of approximately 2fco, which originate from the low 
energy ring structure and, thus, are analogous to the 
Friedel oscillations in fermion systems. 

The spin density, defined as S(r, (f>) = 
V'a( r ) 4>)&aP' l t'p{ r i 4>)i exhibits an interesting topo- 
logical spin texture configuration. Let us first look at its 
distribution along the a;-axis where the supercurrent is 
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along the y-direction and the spin lies in the xz-plane. 
Explicitly, we express S z {r,<j>) = \{\f{r)\ 2 — \g{r)\ 2 ) and 
S x (r,<p) = f(r)g(r). The radial oscillations of |/(r)| 2 
and |g(r)| 2 have an approximate 7r phase shift, which 
arises from the different angular symmetries. As a 
result, S spirals along the x-axis as plotted in Fig. [2] B 
at the pitch value of the density oscillations. The spin 
density distribution in the whole space can be obtained 
through a rotation around the z-axis, which exhibits the 
skyrmion configuration. 

The recent research focus of the "synthetic gauge 
fields" in cold atom systems provides a promising method 
to observe the above proposed exotic spontaneous half- 
quantum vortex [10h12[ 1 141 ]. Artificial magnetic fields 
have been synthesized by light-atom interaction to gen- 
erate vortices experimentally [IH • Similar methods to 
generate SO coupling have been proposed [Io|, 13 , 31 1 , 



by the matrix form as 



and are currently under experimental implementation. A 
convenient method is the "tripod scheme" : three atomic 
degenerate hyperfine ground states are coupled to an ex- 
cited state by laser fields to create a pair of degener- 
ate dark states denoted with pseudospin indices | f) and 
| I) . The spatially dependent laser fields couple the pseu- 
dospin and the momentum of the atoms. A typical set of 
parameter values are provided in Ref. for a trap fre- 
quency lot = 2TrxlOHz,a characteristic Rabi frequency 
= 2u x 10 7 Hz, and a detuning A — 2nx 10 11 Hz, which 
satisfy lot <C £1 <C A and lot <C £1 2 /A. The value of a 
can vary from up to the order of several tens. With the 
typical particle number Nq = 10 6 , the interaction energy 
scale E int is around lOOnK, or 2kHZ [32j, and thus j3 is 
of the order of several tens. Hence, the parameters used 
in the Fig. [5] A and B are realistic in real experiments. 
The experimental detection would be straightforward as 
performed in previous time-of-flight imaging in vortex 



experiments }33l |34| . By separately imaging the density 
profiles of the pseudospin up and down components, the 
radial oscillation of the S z component can be directly 
seen. 

Another class of boson systems exhibiting SO coupling 
is the indirect excitons in 2D coupled double quantum 
well systems. We only need to consider the heavy hole 
(hh) band with the effective mass rn* hh and j z = ±|fi., 
which is separated from the light hole band with a gap of 
the order of lOmeV. The Hamiltonians of the electron- 



hole system reads H e 



Ira 



- V 2 +ih\ R (d etX (Jy - d e ,y<J: 



Hhh — —- 



-Vl, and H e 



m* is the effective mass of the conduction electrons; \r is 
the Rashba SO coupling strength of the conduction elec- 
tron; e is the dielectric constant; and d is the distance 
between the double wells. The center-of-mass motion 
of excitons in the BEC limit can be separated from the 
electron- hole relative motion. Similarly to Ref. (35|, the 
effective Hamiltonian of the 4-component hh excitons de- 
noted as (s e ,jhh) = (±^, |), (±5, —|) can be represented 



where 



H„ 



( E ex (k) H so (k) 

H* ao (k) E ex {k) + A{k) W(k) 

W*{k) E ex (k) + A(k) H ao (k) 

V H* so {k) E ex (k) ) 



'h 1 1 



where k is the center-of-mass momentum; M = m* 
is the total mass of the exciton; E ex (k) = h 2 k 2 /(2M); 
H so {k) = — ^Tj-Xjt(k y +ik x ); A(k) is the exchange integral 
and W{k) has the d-wave structure as (k x + iky) 2 , both 
of which are exponentially suppressed by the tunneling 
barrier between electron and hole layers, and can be ne- 
glected. Consequently, H ex becomes block-diagonalized. 
We consider using circularly polarized light to pump the 
exciton of (— 5, |), and then focus on the left- up block 
with heavy hole spin j z = |. We will then end up with 
the 2D version of the Hamiltonian given by Eq. [TJ 

In summary, we find that bosons with SO coupling ex- 
hibit a novel type of BEC beyond Feynman's argument. 
In the limit of large SO coupling, the degeneracy of the 
single particle ground states in momentum space leads to 
two fragmented condensations with equal particle num- 
ber partition in homogeneous space. In the presence of a 
harmonic trap, the condensate spontaneously breaks TR 
symmetry and develops the half-quantum vortex config- 
uration with the skyrmion type topological spin density 
configurations. 

C. W. is supported by NSF under No. DMR- 0804775, 
and ARO-W911NF0810291. 
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